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ABSTRACT 



We have used the Very Long Basehne Array to measure the trigonometric parallax of several member stars of the Orion Nebula Cluster 
showing non-thermal radio emission. We have determined the distance to the cluster to be 414 ± 7 pc. Our distance determination 
allows for an improved calibration of luminosities and ages of young stars. We have also measured the proper motions of four cluster 
stars which, when accurate radial velocities are measured, will put strong constraints on the origin of the cluster. 
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1. Introduction 

The rich Orion Nebula Cluster, the best-studied of all star clus- 
, ters, is an important "laboratory" for stellar astrophysics and 

■ early stellar evolution. Most of its x 3500 stars formed within the 
_ past 2 10*' yr The most massive stars, the 0' "Tr apezium" sys- 
. tern, have evol ved close to the main sequence dHillenbrandl 1 9971: 

■ iPaUa & Stabler 19 99). while many lower mass members are still 
I pre-main sequence (PMS) stars. The dominant Trapezium star 

0'C (spectral type 05-07) is responsible for the excitation of 
the spectacular Orion Nebula (ON) Messier 42, the only ionized 
nebula visible with the n aked eye and e xtensively studied since 
the seventeenth century ( lHerczeglll998h . since the 1950s at all 
wavelengths. 

The youngest stars in the region are found in the 
nearby Becklin-Neugebauer/Kleinmann-Low (BN/KL) region. 
The BN/KL region is a very dense part of Orion Molecular 
' Cloud - 1, which is located j ust a fraction of a pars ec behind the 
; ONC (IZuckerman 1973; Genzel & Stutzkilll989l) and is home 
, to at least one high-mass st a r that is currently form i ng (Orion- 
l) (IChurchweUetal.1 Il987t iMenten & ReidI 119951: iReid et al.l 

Given the importance of the ONC and the BN/KL region for 
the study of star formation and early stellar evolution, knowl- 
edge of their distance, D, is of great interest. For example, lu- 
minosities are a critical constraint on stellar evolution models, 
which deliver stellar masses and cluster ages. Note that a lumi- 
nosity determination depends on the square of the distance! In 
addition, accurate distances are needed to determine binary-star 
separations, to convert angular motions to velocities, and to esti- 
mate stellar mass-loss rates. 

Distance estimates for the ONC from the early 1900s ranged 
from 2000 pc based on apparent magnitudes and colors of faint 
stars in the ON's neighborhood, to 180 pc, based on the mov- 
ing cluster method using stars wide ly distributed around the 
ON (iPickerind [19171: lKaptevnlll918h . Later, luminosity -based 
distance estimates of the Trapezium stars were plagued by the 
poorly known and spatially variable extinction within the nebula 
and ranged from 540 pc to 300 pc (iTrumplerll 19311; iMinkowskil 



|1946l) . Various distanc e estima t es fro m the 1960s through the 
1980s, summarized in iJeffriesI (l2007h . used optical and near- 
infrared (NIR) photometry and color-magnitude diagrams to 
yield 347 pc < D < 483 pc. 

Uncertain extinction corrections and variations of extinction 
across the Nebula, quantified, e.g., by Bohlin & Savag e ( 198l|), 
may have led to part of the large spread in distances pub- 
lished over the last few decades. Thus, measuring distance by 
methods that are independent of photometry is highly desirable. 
Comparing radial velocities and proper motions for samples of 
stars for which both cou ld be measured has yielded distances 
from 380 pc to 520 pc (|Johnsonlll965t IStrandlll958l) . Recent 
modeling of the rotational properties of ONC pre-main-sequence 
stars has given D = 440 + 34 pc for a sample of 74 stars, which 
falls to 392 + 3 2 pc when stars with accretion disks are excluded 
(|Jeffriesll2007b . Very Long Baseline Interferometry (VLBI) de- 
terminations of the internal proper motions of "clouds" of wa- 
ter vapor (H2O) masers in the BN/KL region led to another 
photometry-indepe ndent distance estimate of D = 480 + 80 pc 
(iGenzel et al.llT98TI) . 

Very recently, iKraus et all (|2007|) used multiple visual and 
NIR interferometric observations to determine two orbital solu- 
tions for the 0'C close binary system. One invokes a distance of 
434 + 12 pc and the other 387 + 1 1 pc. The existing data do not 
allow favoring one over the other. 

It is apparent from the above summary that there is disagree- 
ment of at least +10% among the various methods of estimating 
the distance to the ONC, which leads directly to uncertainties of 
+20% for luminosities. Clearly a much more accurate distance 
for the ONC would be a fundamental advance for stellar astro- 
physics. 

The "gold standard" for astronomical distance measure- 
ments is the trigonometric parallax, which uses the ancient sur- 
veying technique of triangulation. For astronomical applications, 
the Earth's orbit is used as the length scale for one leg of the 
triangle, and the angular displacement of an object (relative 
to very distant sources) directly yields the object's distance. 
Optically, only a single ONC member star has a trigonometric 
parallax, 361^g^'^ pc obtained from the Hipparcos astrometric 
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Star 


'3'J2000 


<^J2000 


STMRAD 


H97 JW88 


COUP 


GMR 














A 


05 35 11.80318 


-5 21 49.2504 








450 


12 


05 35 15.82615 


-5 23 14.1296 


e' AiOri 


680 




745 


G 


05 35 17.95028 


-5 22 45.5058 


MT Ori 


815 


567 


932 


F 


05 35 18.37143 


-5 22 37.4342 


VI 229 Ori 


852 


589 


965 



Table 1. Non- thermally emitting stars in the ONC detected 
by the VLBA. The first c olumn gives the radio source des- 
ignation from iGaray et"an (Il98 7) and the second and third 
columns the J2000 position determined by us and calculated for 
the midpoint of our observations (epoch 2006.46). The iden- 
tifications in the fourth through seventh columns come from 
the SIM BAD database and Hill enbrand (1997.) : , Jones & Walker, 
(ll988h : lGetman et al.l(l2005h . 



satellite (iBertout et al.lll999h : but this measurement is not accu- 
rate enough for most astrophysical applications. Extraordinarily 
high astrometric accuracies can now be achieved with VLBI 
techniques. Recently, trigonometric parallax distances to the ON 
have been reported to be 437 + 19 pc from water masers in the 
BN/KL region usin g the VERA (VL BI Exploration of Radio 
Astrometry) array (iHirota et al.l l2007l) and 389^2i_P'^ from an 
ONC star (GMR A; see below) using the NRACQ Very Long 
Baseline Array (VLBA) dSandstrom et alJ l2007l) . These mea- 
surements represent a significant advance for the distance to the 
Nebula, but the values still differ by 48 + 30 pc. 

Using the VLBA, source positions can be now be measured 
with an accuracy of ~ 10 micro-arcseconds (/ias) relative to very 
distant quasars. Recently, VLBA observations of methanol and 
water masers have yielded trigonometric parallax distances for 
the star-forming region W30H in the Perseus s piral arm of the 
Milky Way of 2.0 kp c with w 2% accuracy (iXu et al.1 l2006t 
iHachisuka et al.l 12006*). After completing these measurements, 
we started an astrometric program with the VLBA in order to 
determine trigonometric parallaxes for radio-emitting stars in the 
ONC. 

The Orion Nebula is known to contain about 90 com- 
pact radio stars, either associated with the ONC o r with dust- 
embedded sources in the ne arby BN/KL region i Garav et~an 
ll987HChurchwell et al.ll987l) . The emission from many of these 
sources is thermal in nature, either resulting from external ion- 
ization of circumstellar material (protoplanetary disks or "pro- 
plyds") by the UV radiation of fl'C Ori, or from internal UV 
sources exciting hypercompact ionized regions. However, more 
than 30 sou rces have shown significant time v ariability and/or 
emit x-rays (" Felfi et al.ll99llZapata et alJ2004 ). Probably these 
are non-thermal radio emitters, making them suitable sources for 
VLBA parallax measurement^ 

2. Observations and data analysis 

2.1. Observations 

Our observations with the VLBA were conducted on 2005 
September 25, 2006 March 02, 2006 September 09, and 2007 

' The National Radio Astronomy Observatory (NRAO) is operated 
by Associated Universities, Inc., under a cooperative agreement with 
the National Science Foundation. 

^ The VLBA images containing stars presented in this paper have 
an rms noise level around 0.1 miy beam"' . In a 2 x 1 milliarcsecond 
FWHM synthesized beam this corresponds to a brightness temperature 
of 8.6 10' K. 



March 6. These dates well sample the maximum extent of the 
Earth's orbit as viewed by the source. We designed the obser- 
vations to measure only the Right Ascension (R.A.) component 
of the trigonometric parallax signature, because this component 
has about twice the amplitude of the Declination (Decl.) compo- 
nent and the angular resolution of the VLBA for this source is 
about two times better in R.A. than in Decl. 

The pointing position for the VLBA antennas for the Orion 
Nebular Cluster was (a, 5)j2ooo = 05'^35'"15."0, -05°22'45". 
We chose the compact extragalactic source J054 1-054 
as the phase-reference from the Internat i onal Celestial 
Reference Frame (ICRF) catalog dMa et all Il998l) . adopt- 
ing its position to be {a,5)j2ooo = 05Ml"38."083384 + 
0.'000019,-05°41'49:'42839 + 0."00046. We alternately ob- 
served J054 1-054 and the Orion position, switching sources 
every 40 seconds over a period of 8 hours during each epoch. 
The observing setup used 8 intermediate frequency (IF) bands 
of 8 MHz, with 4 detecting right and 4 left circular polarization. 
The center frequency of the bands was 8437 MHz and togther 
they spanned a total 32 MHz in each polarization. 

The data were correlated with the VLBA correlator in 
Socorro, NM, with an averaging time of 0.131 seconds. The 
position offsets of the ONC stars from the pointing position are 
small relative to the FWHM of an individual VLBA antenna's 
primary beam (5.'3). However, we had to correlate the data 
in three passes, owing to the large interferometer fringe-rate 
differences among them. Three correlator phase-center positions 
were used (a,6)moo = 05''35"1 1.'8022, -05°21'49:'229, 
05''35'"16.'2890, -05°23'16:'575, and 05*'35'"18.'3706, 
-05°22'37!'436. Bandwidth smearing limits the effective field- 
of-view to about 22 arcseconds about each coiTelator position. 
Within this field of view, amplitude decoiTelation is less than 
10% on our longest interferometer baselines. 

2.2. Calibration 

Most of the data calibration steps are standard and described in 
the on-line documentation of the NRAO AIPS software pack- 
ag41 Below we discuss additional calibrations that allow im- 
proved astrometric accuracy. 

The main source of systematic error for cm-wave phase- 
referenced observations is uncompensated interferometric de- 
lays introduced by the Earth's atmosphere and ionosphere, and 
we conducted supplem entary observations that allowed us to 
correct for these effects jReid & Brunthalerll2004l) . These obser- 
vations were done by spreading the 8 left-circularly polarized IF 
bands to sample about 500 MHz of bandwidth. Approximately 
15 ICRF sources whose positions are known to better than 1 mas 
were observed over a span of about 40 minutes at the beginning, 
middle and end of the Orion observations. The effects of iono- 
spheric delays were removed by using total electron content val- 
ues from GPS data. Also, data were corrected for the best Earth's 
orientation parameter values, as the VLBA correlator by neces- 
sity uses preliminary values available at correlation. Broad-band 
interferometer residual delays were estimated and modeled as 
owing to clock drifts and zenith atmospheric delays. The phase- 
referenced data was then coiTected for these effects. 

Electronic delay and phase differences among the different 
IF bands were removed by measuring them on a strong cali- 
brator, 0530H-133, and coiTecting the data. Then the data on the 
rapid-switching reference source, J0541-054, from all IF bands 
were combined and tables of antenna-based interferometer phase 

^ See the AIPS "cookbook" under http://www.nrao.edu/aipsl 
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versus time were generated. These phases were interpolated to 
the times of the Orion data and subtracted from that data. 

All calibrations were applied to the data and synthesis maps 
were made using the AIPS task IMAGR. Typical dirty-beam 
FWHM sizes were 2.1 by 0.9 mas elongated in the north-south 
(NS) direction. Positions of both the reference source and the 
Orion stars were determined by fitting 2-dimensional Gaussian 
brightness distributions to the images using JMFIT. 

2.3. Imaging 

We searched for compact radio emission from a total of 23 
stars, but only detected the four stars in Table 1 above a se- 
cure (> 5cr) detection threshold of w 1 mJy. The locations of 
these stars within the core of the ONC are shown in Fig. 1, us- 
ing the stellar des ignations from Garay, Moran & Reid (GMR, 
iGarav etaTlllQSTh . The star GMR A has no optical counterpart, 
but has long been known to be a variable radio source, whose 
emission was found by VLBA observations to be very compact: 
~ 1 milli-arcsecond (mas). In 2003 January a remarkable out- 
burst of several days duration was observed at millimeter and 
radio wavelengths, during w hich the source's 8 6 GHz flux den- 
sity rose more than 10-fold dBower et al.l 12003 ). Near infrared 
(NIR) photometry and spectroscopy indicate that GMR A is a 
deeply embedded, very young (T Tauri) star with more than 20 
mag of visual extinction. In x-rays, GMR A is characterized as 
a "flare source" and the x-ray data indicate an absorbing column 
of 4 10^^ cm"^, w hich is consistent with the extinction derived 
from the NIR data (iFeigelson et al.ll2002h . 

The other stars detected, GMR 12, GMR F and GMR G, 
are coincident with optically visible stars and are undoubtedly 
ONC members. The later two stars are listed as ;r l925 and 
7rl910, respectively, in the classic Parenago catalog ("Parenagd 
[19541) . GMR 12 is the binary companion of the Trapezium star 
0' A Ori (named 0' A2 Ori) . This 0."2 separation binary sys- 
tem has recently been resolved by N IR speckle interferometry 
jPetr et al.lll998i: IWeigelt et al.lll999l) . (The primary, 0^ Aj Ori 
(BO. 5 V), itself is a spectroscopic binary with a cooler compan- 
ion.) GMR 12 has been detected previously wi th VLBI observa- 
tions jFelli et al.lll9"89ill99lHGaiTington et al.l l2002). 

Three of our four stars were detected at all four epochs; 
GMR G was not detected on the first epoch owing to source vari- 
ability. Our observations delivered the positions listed in Table 
1 with an estimated precision of about 0.13 mas in R.A. and 
0.28 mas in Decl. relative to the extragalactic reference source. 
Absolute position uncertainties are dominated by the ^ 0.5 mas 
uncertainty in the ICRF position of the reference source. Our 
position s are in excellent agreement with the values determined 
by Gomez et al.l (12005) using the Very Large Array within their 
larger errors. 

Single component JMFIT solutions of our reference source, 
0541-054, for each epoch yield an elongated elliptical Gaussian 
component with position angles, PAs, between 145 and 158 de- 
grees (east of north) that is barely resolved. JMFIT delivers de- 
convolved sizes between 0.8 and 1 .6 mas for the major axis and 
between 0.3 and 0.4 for the minor axis. This evokes the possi- 
bility of a core-jet structure. Note, that to first order, an evolving 
jet would contribute to the proper motion, but not to the paral- 
lax. Further note that our parallax essentially derives from the 
east-west (EW) data. 

To quantitatively study the effect that source structure might 
have on the determinations of 0541 -041 's centroid position we 
conducted 2 (Gaussian) component fits using JMFIT. In this case 
the task delivers a very compact stronger component of size 



< 0.7 mas and a weaker elongated component with a major axis 
of size a! 4 mas. For the latter' s minor axis JMFIT returns up- 
per limits between 0.05 and 0.2 mas for the first three epochs, 
while it has broadened to 2 mas in the forth epoch. This "broad- 
ening" may indicate more complex source structure emerging. 
For the compact component, which almost certainly represents 
the QSO's core, we found small deviations from the 1 Gaussian 
fit centroid position of -0.057, -0.037, -0.064, and -1-0.005 mas 
in R.A. direction and -0.084, -1-0.036, -1-0.002, and +0.002 mas 
in Decl. direction for the four epochs, which contribute to the 
error floors discussed in ^ 

In Table 2, we present the stars' positions together with the 
fitted peak intensities, Sp, (brightness values) and integrated flux 
densities, S . A comparison of the numerical values of 5p and S 
shows that all sources at most epochs appear only marginally re- 
solved, if at all. Since the sources appear more resolved when 
they have low flux densities and do hardly appear resolved at all 
when their flux densities are highest, we ascribe this apparent re- 
solving of sources to residual calibration uncertainties that affect 
low signal-to-noise cases relatively the most. 



3. Results - Parallax and proper motion 
determinations 

We used our positions as input data for a least-squares fitting 
program, which modeled the data as the sum of the sinusoidal 
parallax term and a proper motion term. The parallax term is en- 
tirely determined by one parameter (tt), since the Earth's orbit 
and source directions are well known. The proper motion term 
requires two parameters for each coordinate: an angular offset 
and speed. The first four lines of Table 3 give the best fit paral- 
laxes and proper motions for all four stars individually. In Fig. 
2 we plot the positions of the 4 stars versus time with the best 
fit parallax and proper motion models. Since all stars are at the 
same distance, within measurement accuracy, we also present a 
combined solution in which we simultaneously solve for a sin- 
gle parallax parameter and separate proper motions for the three 
stars detected in all four epochs (see final line of Table 3 and Fig. 
3). 

For this fit, an error "floor" was added to the formal position 
uncertainties in order to account for systematic errors, resulting 
mostly from unmodeled atmospheric delays and variations in the 
calibrator positions; see 32.21 The error floors for the R.A. and 
Decl. data (0.12 and 0.27 mas, respectively) were separately ad- 
justed to bring the reduced value (per degree of freedom) to 
unity in each coordinate. Without the error floors the reduced 
values were 3.7 in R.A. and 7.5 in Decl. (per degree of freedom), 
respectively. After applying these weighting factors, the best fit 
parallax was 2.415 ± 0.040 mas, corresponding to a distance of 
414.0 ± 6.8 pc. Both approaches, the single source average and 
the joint solution fit, give similar results, and we adopt the joint 
solution value, i.e., a distance of 414 ± 7 pc. 

We also determined a parallax including the data for GMR G. 
For this fit, we obtained n = 2.448 + 0.042 mas, corresponding 
to D = 408.5 + 7.0 pc. This is within Icr with the 414.0 pc 
determined above, a value which we retain as our best estimate 
for the reason that of all sources GMR G consistently showed the 
lowest flux densities making it the most vulnerable to systematic 
errors. This is reflected by the fact that the formal error of the fit 
that included GMR G is not smaller than that of the three source 
only fit. 

Our distance, 414 pc, is the most accurate measurement for 
the Orion region, having an uncertainty of only 1 .7%. Both of the 
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Star 


Epoch 


Sp 




5 




Q'J2000 










PA 






(mJyb-') 




(mJy) 














(mas) 


(mas) 


(deg) 


GMR A 


1 


0.9(0.1) 


0, 


.9(0.3) 


05 


35 


11 


.803276(6) 


-05 21 49, 


.24884(20) 


< 1.6 


< 0.8 


161(3) 




2 


2.1(0.1) 


2, 


,8(0.3) 


05 


35 


11 


.803031(2) 


-05 21 49, 


.25059(6) 


1.1(3) 


0.6(2) 


171(30) 




3 


10.7(0.1) 


13, 


,1(0.2) 


05 


35 


11, 


,803409(1) 


-05 21 49, 


,25063(2) 


1.5(2) 


0.5(1) 


170(4) 




4 


10.8(0.1) 


14. 


,9(0.2) 


05 


35 


11, 


,803135(1) 


-05 21 49, 


,25219(1) 


0.8(1) 


0.4(2) 


71(8) 


GMR 12 


1 


2.1(0.1) 


6 


,9(0.5) 


05 


35 


15, 


,826115(5) 


-05 23 14, 


,12875(9) 


3.1(3) 




58(7) 




2 


14.4(0.1) 


16, 


,8(0.2) 


05 


35 


15, 


,825945(1) 


-05 23 14, 


,12943(1) 


1.0(1) 


0.3(1) 


163(3) 




3 


4.9(0.1) 


6, 


,8(0.2) 


05 


35 


15, 


,826430(1) 


-05 23 14, 


,12961(3) 


1.6(2) 


0.7(1) 


173(6) 




4 


4.8(0.1) 


5, 


,5(0.2) 


05 


35 


15, 


,826262(1) 


-05 23 14, 


,13120(2) 








GMRG 


1 


< 0.24 


























2 


1.1(0.1) 


1, 


,5(0.2) 


05 


35 


17, 


,950071(2) 


-05 22 45, 


,45702(9) 


1.4(4) 


< 0.7 


19(20) 




3 


2.0(0.1) 


2, 


,2(0.2) 


05 


35 


17, 


,950564(2) 


-05 22 45, 


,45477(7) 


< 1.5 




9(40) 




4 


4.7(0.1) 


5, 


,1(0.2) 


05 


35 


17, 


,950357(1) 


-05 22 45, 


,45364(2) 


< 0.7 


< 0.6 


34(25) 


GMRF 


1 


1.0(0.1) 


1, 


,9(0.2) 


05 


35 


18, 


,371521(4) 


-05 22 37, 


,43477(14) 


2.1(5) 


< 0.9 


162(9) 




2 


6.8(0.1) 


7, 


,3(0.2) 


05 


35 


18, 


,371255(1) 


-05 22 37, 


,43470(2) 


< 0.6 


< 0.4 


16(43) 




3 


26.8(0.2) 


29, 


,2(0.3) 


05 


35 


18, 


,371659(1) 


-05 22 37, 


,43379(1) 


1.4(1) 


< 0.1 


167(3) 




4 


2.4(0.1) 


3, 


,1(0.2) 


05 


35 


18, 


,371415(1) 


-05 22 37, 


,43403(3) 


< 1.0 


< 0.9 


124(35) 



Table 2. Fitted values for the four ONC stars determined by Gaussian fitting. Left to right we list name, epoch, peak intensity 
(brightness), total flux density, R.A., Decl., fitted major and minor axis and position angle. Formal Icr uncertainties in the last listed 
digit returned by JMFIT are given in parentheses. For R.A. values with smaller uncertainties than 1 //second, 1 is given for display 
reasons. Epochs 1 through 4 correspond to 2005 September 25, 2006 March 2, 2006 September 9, and 2007 March 6. Absolute 
position accuracy is limited the +0.5 mas ICRF position accuracy of the phase-reference source, J0541-054. GMR G was not 
detected at epoch 1 . The quoted brightness is a 3cr upper limit. For some fits JMFIT did not return a value for the major or minor 
axis, for some only an upper limit. 



trigonometric parallax es mentioned above (iHirota et al.l 120071 
ISandstrom et al.ll2007l) are consistent with our value within their 
larger uncertainties of about +20 pc. 

4. Discussion 

Our accurate distance is approximately 10% lower than the 
450 pc distance often assumed for the ONC. Since luminosity is 
proportional to distance squared, our distance implies that lumi- 
nosities of the ONC member stars have been over-estimated by 
nearly 20%. Reducing luminosities affects theoretical investiga - 
tions of the ONC's star formation history ( Palla & Stahled 19991) . 
Age and mass estimates from stellar evolution mod els are dis- 
cussed in the literature. We have used an online tool (ISiess et al.l 
|2000; Siess 2001) which allows one to calculate the change in 
age and mass, M, for pre-main-sequence stars with masses be- 
tween 0.1 and 7 Mq and to trace lines for given masses in a lu- 
minosity versus temperature diagram ("PMS tracks"). For low- 
mass (M < 1 Mq), cool (< 5000 K) stars, the vertical tracks 
indicate little change in mass from a reduced luminosity and we 
find little (at most a few percent) difference in mass for a given 
luminosity. However, ages are significantly affected and become 
20 to 30% larger by using our accurate distance. The same hap- 
pens for higher mass stars, except that masses are affected, too, 
and drop by up to 10% for stars with 7 Mq. 

The proper motions of the four stars are in the heliocentric 
equatorial reference frame. We transform them into a Galacto- 
centric Cartesian reference frame (U, V, W), where U points 
toward the Galactic center, V points toward the direction of 
Galactic rotation, and W points toward the North Galactic Pole. 
To correct for the Sun's motion, we us ed the latest values derive d 
from the analysis of Hipparcos data dPehnen & Binnevlll998l) . 
For the distance to the Galactic center and the circular rotation 
speed in the Solar Neighborhood, we assume the lAU values of 
8.5 kpc and 220 km s"', respectively. 



Source 


n 


D 








(mas"') 


(pc) 


(mas y"') 


(mas y ' ) 


GMR A 


2.390(0.104) 


418.4(18.2) 


+ 1.82(0.09) 


-2.05(0.18) 


GMR 12 


2.393(0.053) 


417.9(9.2) 


+4.82(0.09) 


-1.54(0.18) 


GMRG 






+4.29(0.17) 


+3.33(0.37) 


GMRF 


2.462(0.051) 


406.1(8.4) 


+2.24(0.09) 


+0.66(0.18) 


Mean 


2.425(0.035) 


412.4(6.0) 






Joint 


2.415(0.040) 


414.0(6.8) 







Table 3. Fitted distances and proper motions. The first col- 
umn lists the radio source designation. The second and the third 
columns give the parallax and distance estimates derived from 
a fit to the measured positions in which the parallax and the 
proper motions in R.A. and Decl. direction (listed in columns 
four and five) were left as free parameters. A proper motion of 
1.0 mas y"' corresponds to a velocity of 2.0 km s"' at a distance 
of 414 pc. The fifth line gives the weighted mean parallax and 
distance and their standard errors of the mean, and the sixth line 
gives the results of a combined solution, where only one parallax 
was fitted to the measurements of GMR A, GMR 12, and GMR 
F. 



The calculated U, V, and W values are dependent on the as- 
sumed radial velocities of our stars. Unfortunately, no accurate 
radial velocity values have been published for GMR 12, F, and G. 
With the helioc entric radial veloci ty determined for GMR A of 
+ 14 + 5kms"' dBower et alj2003h . we calculate U = +9.6+4.2 
km s"'. Ay = -3.8 + 2.6 km s"' (slower than Galactic rotation), 
and W - +5.5 + 1.7 km s"', where AV is the difference between 
V and the rotation speed of the Milky Way at the position of the 
ONC. For a rotation curve that is constant with Galacto-centric 
radius, as appears to be the case near the Sun, the assumed value 
of the rotation speed does not significantly affect AV. 
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If accurate radial velocities for all our program stars could 
be determined, e.g., with IR spectroscopy, such measurements 
combined with our proper motions could constrain dynamical 
scenarios for the formation of the Orion complex, which is lo- 
cated at the large distance of 19.4 degrees (140 pc) below the 
Galactic plane. To explain this puzzling offset one theory invokes 
a high velocity cloud colliding with the Galactic plane 6 10^ y 
ago, sweeping up and compressing mate rial in the process, and 
subsequently oscillating about the plane (iFranco et alJll988h . 

The proper motions of our four stars have a dispersion of 1 .5 
and 2.6 mas y ' in R.A. and Decl. These values, which corre- 
spond to 3 and 5 km s"', respectively, are somewhat uncertain 
because of small number statistics, but are significantly larger 
than velocity di spersions of 1 km s~' found from optical prope r 
motion studies ("J ones & Walked 1 19881 1 van Altena et al. I [Mi- 
One explanation for this difference may be that three or our four 
stars lie far from the cluster center and, thus, could be respond- 
ing to a larger enclosed mass than the bulk of the stars nearer to 
the center. 
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Fig. 1. Infrared K-band image of the ONC and the BN/KL region, (courtesy M. McCaughrean). The Trapezium stars (A, B, C, D & 
E) and the BN object are labeled. Our program stars (OMR A, OMR F, OMR G & GMR 12) are labeled and encircled. The image 
was taken with the ESO Very Large Telescope. 
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Fig. 2. Position versus time for GMR G (black), GMR F (red), GMR 12(green), and GMR A (blue). The top and bottom panels 
show the eastward (R. A. cos(Decl.)) and northward (Decl.) offsets, respectively. For each star the best fit distance from the correlator 
phase-center position has been removed and then offset for clarity. Also plotted are best fit models. For each star, except GMR G, 
the model is characterized by five parameters: one for the parallax and two for the proper motion speed and position offset for each 
coordinate. For GMR G, which was detected at only 3 epochs, the parallax value was fixed at the "joint- solution" value (2.415 mas; 
see Fig. 3) and only the proper motion parameters were adjusted. 
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Fig. 3. Joint parallax fit for the 3 stars detected at all 4 epochs using position versus time data for GMR F (red), GMR 12(green), and 
GMR A (blue). The top and bottom panels show the eastward and northward offsets, respectively. For each star the best fit proper 
motions have been subtracted. The best fit parallax was 2.415 + 0.040 mas, corresponding to a distance of 414.0 + 6.8 pc. 



